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ABSTRACT 

Fine-structure mid-infrared emission lines with critical densities in the regime 10 4 ^ 
n crit ^ 10 6 cm -3 can be employed to probe the outflow from Wolf-Rayet stars at radii 
of ~ 10 15 cm. Narrow-band mid-infrared imaging and spectroscopy of the nearest WR 
star to the sun, 7 Velorum is analysed for spatially resolved forbidden line emission in 
the Wolf-Rayet outer wind. The [Siv] 10.52-/mi and [Neil] 12.81-/xm emission regions 
are found to be spatially extended, compared to unresolved continuum and He and 
C recombination line emission. The [S iv] and [Ne 11] emission line distributions have 
a high degree of azimuthal symmetry, indicating a spherically symmetric outflow. A 
model wind with a modest degree of clumping (clumping factor / ~ 10) provides a 
better match to the observations than an unclumped model, The overall line inten- 
sity distributions are consistent with a freely expanding, spherically symmetric 1/r 2 
outflow with constant ionization fraction and modestly clumped density structure. 

Key words: circumstellar matter - infrared: stars - stars: Wolf-Rayet - stars:winds, 
outflows - stars: individual: 7 Velorum 



1 INTRODUCTION 

7 Velorum is a binary sys tem comprising a WC8- 
type Wolf-Rayet JSmith 1968i ^ and 07.5 III companion 



jPe Marco fc Schmutj 1 19991 '). The bi nary is well-defined , 



with a period of 78.53 days, i ~ 65° ( Sch mutz et al 
and a separation of ~3.8 milliarcsec (|North et al 



1997) 



2007 



iMillour et ai1l2007l ). The O-star appears typical and is the 
more massive st ar in the binary with Mp = 28.5Mp) against 
M« B = 9.0M^ dConti fc Smithll 19721 ; |Pe Marco fe Schmutj 



1 19991 ; iNorth et all 120071 '). TheWC8 is the nearest example 
of a Wolf-Rayet st ar, at a distance of 342 pc (see below; 
Ivan Leeuwenll2007l ). affording the greatest opportunity to 
study this important class in depth 

The outflow from the WC8 is considerable, with a mass- 
loss rate of 2.95xl0" 5 sC M yr _1 sC 8.85xl0" 5 if un- 
clumped (and a factor of 3.8 lower if clump ed by a factor of 
10), at least 100 times that of the O-star ( |Pe Marco et al.l 
l200d ; Barlow, Roche & Aitken 1988). The ejecta can be de- 
tected at great separations from the star, several times the 
binary separation, from the Brem sstrahlung continuum at 
radio wavelengths JSeaa,uistlll976l) and from l ow-density for- 
bidden emission lines (e.g. Barlow et all 198*81 ). The wind ter- 
minal velocity, found from those forbidden emission lines, is 
~1520 kms -1 , confirmed by UV eclipse absorption line pro- 
files (St-Louis, Willis & Stevens 1993). An energetic wind- 
collision zone is detected in UV and X-rays, with variability 
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consistent with t he effect of the alignment of the binary 
to the observer (|St-Louis et all Il993l ; IStevens et ail 11999 ; 
ISchild et aj|2004l ). X-ray absorption column densities sug- 
gest a lower mass-loss rate by an order of magni tude, which 
may be evidence of clumping in the inner wind (|S child et all 
120041) . 

Mid-infrared spectroscopy of 7 Vel has been used to 
investigate the wind structure. One important motivation 
was to determine the neon abundance in the wind, as Wolf- 
Rayet stars are thought to be the origin of most 22 Ne in 
the ISM. This was done from the forbidden line emission 
[Ne ii ] at 12.81-/im and [Nem] at 15.56-/xm l|Barlow et al.l 
1988, van der Hucht et al 1996) both of which have wide, 
flat spectral profiles suggesting that much of the emission 
arises from the outer regions where terminal velocity has 
been reached. However, from an uncertainty in the mass-loss 
rate due to the un certain distance, th e Ne enrichment was 
not initially found (Barl ow" et al.l fl988). Further analysis in 
the mid-IR and the introduction of wind-clu mping models 
suggests that Ne is enriched in WC star winds ( Dessart et al.l 
|2000| ; ISmith fc Houck1l2005l ; llgnace et al.ll20071 ). 

In order to correctly calibrate the spatial scale of the 
observed outflow it is necessary to fix the distance to the 
object. The distance to 7 Vel has been the subject of 
many revisions in recent years. In this paper we use the re- 
vised Hipparcos parallax of 2.92±0.3 mas from Ivan Leeuwenl 
(2007). This translates to a distance of 342^31'g pc, which is 
broadly consistent with the interfer ometry-derived dist ances 
of 336 pc and 368 pc estimated bv lNorthet all (|2007l ) and 
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Date 


Pilter 


\ 

"centre 


-^50% 


Ori-sourcc 


CUT) 


name 


(jtm) 


(H 


time (min) 


Imaging 










2006 Apr 10 


[Neil] 


12.81 


12.68-12.91 


34.5 


2006 Apr 18 


[Nelljcont 


13.10 


12.95-13.17 


34.5 


2006 Apr 27 


[Siv] 


10.52 


10.37-10.54 


33.4 


2006 May 06 


[Ar in] 


8.99 


8.93-9.06 


33.4 


Spectroscopy, 


slit PA=0° 








2006 May 22 


N broad 


10.36 




16.3 



Table 1. Log of T-ReCS imaging and spectroscopy of 7 Vel 



iMillour et al.l (|2007l ) respectively. This distance is substan- 
tially greater than the initial Hipparcos estimate of 258 pc. 

In this paper, we present spatially resolved 8-13 /xm 
spectroscopy and narrow-band mid-infrared images of the 
7 Vel outflow and investigate the density structure of the 
wind. 



2 DATA 

The imaging and spectroscopic observations of 7 Velorum 
were obtained with the Gemini Thermal-R egion Camera 
Spectrograph (T-ReCS) (|Telesco et al.lll998h in April and 
May 2006 under Gemini programme GS-2006A-Q-64. The 
instrument uses the telescope in the chop-nod mode, with a 
maximum chop-throw of 15", used for all of the observations 
presented here. The fixed pixel size of T-ReCS is 0.09". Im- 
ages of the WR-O binary and a standard star (SAO 219422) 
were taken in four filters covering the [Ne n] and [S iv] emis- 
sion lines and nearby continuum regions (see Table [T]). Spec- 
tra were obtained through a 0.35" slit and the broad- A^ fil- 
ter, centred at 10.36-/im and with 50% transmission between 
7.70-12.97 /im. The spectral resolution was 0.08/xm, sam- 
pled at 0.022 /im/pixel. Three exposures were taken on the 
night, with the total on-source time given in Table [T] Ob- 
serving conditions were good apart from the night of the 
[Arm] imaging observations, which caused some striping in 
the final image, and the spectra were obtained at rather high 
airmass as the target set. A log of the imaging and spectro- 
scopic observations is provided in Table [T] 

Data reduction was completed using the T-ReCS iraf 
pipeline. The chop frames were mean-averaged and the alter- 
nating positions subtracted to remove residual background 
structure. Examination of the spatial-axis distribution of the 
flux suggested there would be no significant off-source emis- 
sion within the chop throw. The spectra were centred and 
aligned along along one pixel-row, though the corrections 
applied were minimal. A wavelength calibration was applied 
using the sky emission line list supplied by Gemini for the N- 
band. The resulting spectra were then averaged together to 
increase the signal-to-noise ratio of the 7 Vel spectrum. Sim- 
ilar processing was applied to the standard star spectrum of 
HD89998. Telluric-correction and flux-calibration were per- 
formed by dividing the averaged 7 Vel by the standard star 
spectrum and multiplying by a flux-calibrated blackbody to 
match the temperature of the star; 4280 K deri ved from 
the s pectral-type of Klin for HD89998 (SIMBAD; [johnsonl 
1966). A 1-D extraction of the spectrum was taken, sum- 
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Figure 1. Spectra of 7 Velorum at the peak mid-IR flux and 
at spatial separations to the North and South. The [S iv] and 
[Nell] emission line components at 10.5 and 12.8 (im increase in 
prominence away from the central position. 



ming the spatial-axis rows, and, separately, the individual 
rows along the spatial-axis of the 2-D spectrum were binned 
into 3-pixel bins, mean-averaging the rows to produce the 
binned spatial spectra. 



3 RESULTS 

3.1 Spectroscopy 

Cuts across the spatial-axis of the spectrum of 7 Vel yield 
FWHM of 0.45" in the continuum emission between 10 and 
12 ^tm, which matches the spectroscopic standard star PSF 
and indicates an unresolved core emission source. There are 
several prominent emission line blends present in the spec- 
trum. These are listed in Table together with a number 
of weaker emission lines and the principal lines contribut- 
ing to each blend. Five of the most prominent lines are at- 
tributed to helium or carbon recombination line blends, but 
the lines at 10.52 and 12.81 fim have contributions from 
the ground state fine-structure S 3+ and Ne + transitions re- 
spectively, blended with helium and /or carbon, see Tabled 
In the WC8 star, hydrogen lines are not expected and are 
not seen in optical or UV spectra (e.g. Willis 1982). Ionized 
carbon is, of course, an important emission component, with 
C ill and C iv recombinat ion lines previously observed i n the 
near- and mid- IR (e.g. lAitken. Roche fc Alien! [1982] ) and 
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^obs 


Identification 


Flux 


H 




(10- 20 Wcm- 2 ) 


8.25 


Civ 19-17, CIII 28-20 


<10 


8.48 


Cm 16-14, Cm 21-17 


<15 


8.67 


Hei 14-8, Cm 25-19, Civ 15-14 


<20 


8.78 


Hei 10-7, Hen 20-14 


<40 


9.72 


Hen 10-9, Hen 13-11, Civ 20-18 


179.8 


9.86 


Cm 13-12, Hei 18-9, Hen 25-16 


63.2 


10.52 


[Siv], Hei 12-8, Hen 21-15 


54.4 




Hen 24-16, Cm 20-17, Civ 16-15 




11.32 


Hei 9-7, Hen 23-16, Hen 18-14 


66.4 




Hen 16-13, Civ 21-19 




12.38 


Hei 7-6, Hei 11-8, Hen 14-12 


155.6 




Hen 22-16 




12.58 


Hei 14-9 


13.4 


12.81 


[Nell], Civ 17-16 


53.4 



Table 2. Wavelengths of prominent emission lines together with 
the most important components in the blended lines. The A and 
fluxes are from Gaussian fits to lines in the centre bin. Errors on 
the fluxes are estimated to be at most 10% of the flux of the most 
prominent lines. 



UV llCastor fe Nussbaumerlll972l ; Ivan der Hucht fe Lamersl 
1973). The blended line fluxes measured in the peak bin are 
listed in Table [2] and the spatially-binned blended line fluxes 
are provided in Table [3] 

The fact that the recombination line blends exhibit 
comparable spatial profiles to the continuum emission, with 
little dependence on wavelength, indicates that seeing con- 
tributes substantially to the profiles and masks the effects 
of diffraction across the 9.7 to 12.4 /im wavelength range; 
over this range, the diffraction-limited image width increases 
from 0.26 to 0.34 arcsec. The continuum and recombina- 
tion line emission is centrally concentrated from an unre- 
solved emission source, arising from the dense, inner regions 
of the wind. The lines at 10.5 and 12.81-/im show a different, 
spatially-extended emission profile from the stellar contin- 
uum and He lines. The blended and deblended (see below) 
emission is plotted in Fig. [2] and Fig. [4] respectively. 

The [S iv] emission is blended with helium and carbon 
emission at 10.5-/xm and the [Nell] emission is blended with 
Civ 17-16 at 12.77^m (Tabled]). The 10.5-^m blend shows a 
FWHM comparable to the unresolved continuum, reflecting 
the dominance of the He/C contribution to that line in the 
inner region. However, beyond ~0.5" the line flux decreases 
more slowly than the continuum, indicating resolved emis- 
sion from [S iv] at large radii in the wind. The 12.81-/im 
emission is less heavily blended. The largest contribution 
is from [Ne n] while the contribution from C iv emission is 
weak, as can be seen from the small flux in the 8.66-/im line 
- the strongest C iv emission line in this spectral range. The 
12.8/im line blend has a FWHM =0.65", significantly larger 
than those of the continua and 10.5 /im emission line blend, 
which have FWHM in the range 0.43 to 0.48 ". This reflects 
the dominant contribution of [Nell] in the 12.8 fim blend 
and indicates an extended distribution. Careful correction 
for He and C lines at 10.5/im is required to extract the [S 
iv] distribution. 

To deblend the lines, the intrat program (supplied by 
Prof P. Storey) was used to compute line ratios between the 
different helium and carbon transitions, allowing for those 
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Separation from peak (arcsec) 

Figure 2. Normalised spatial profiles of emission lines and con- 
tinuum from 7 Vel. The continuum emission at a number of wave- 
lengths between 10 and 12 (im is shown by the thin solid lines and 
has a very similar profile to the those of the recombination lines, 
illustrated by the Hei 7-6 and 9-7 lines and the Hen 10-9 line. 
These are unresolved, but the lines at 10.5 and 12.8/im, which 
contain contributions from [S iv] and [Ne n] have more extended 
profiles. Negative separations are to the South, positive to the 
North. 

components of the 10.5 and 12.81-/im emission blends to be 
calculated and subtracted. The reference lines against which 
ratios of each species would be computed, were Hei 7-6, 
Hen 10-9, Cm 13-12 and Civ 15-14, each the strongest 
line of their species in the 8-13 fim spectral region, in- 
trat calculates the line flux ratios for hydrogenic species 
in specifi ed densities and temper atures using case B recom- 
bination (|Baker fe Menzej 119381 ), which is justified for the 
opaque inner wind where the lines are thought to form (e.g. 
IStevens et al.|[l996j ). The electron temperature and density 
adopted were 30,000 K and 10 10 cm -3 respectively, but the 
dependence of the line ratios on these quantities is fairly 
weak. The contributions from weak lines were estimated by 
generating a synthetic spectrum from the line ratios gen- 
erated from intrat and fitting to the peak flux spectrum, 
where the He and C emission lines are strongest. The contri- 
butions from Hei, Hell, Cm and Civ were adjusted to give 
a good match to the most prominent emission lines. Once 
the ratios of He I, He II, C ill and C iv to each other were de- 
termined, they were applied to the measured fluxes at offset 
positions (see Fig. |3J). The synthetic spectrum was broadly 
consistent with the offset spectra, giving confidence that the 
He and C line contributions to the 10.52 and 12.81-pm lines 
could be estimated. 

Inspection of Fig. [3] shows that the spectral structure 
between 10 and 12 fj,m, which results from high-ra recom- 
bination transitions is well matched by the synthetic spec- 
trum. The region between 8 and 10 (im is noisier because 
of the poorer atmospheric transmission. Note that the dom- 
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Position 
(arcsec) 


\). i A fim 


10.52 


r lux (1U Wcm 
/zm 11.30 /iin 


) 

12.37 [im 


Iz.ol //m 


-1.62 


0.322 


1.023 


0.493 


0.977 


0.680 


-1.35 


1.598 


1.003 


0.567 


0.249 


1.347 


-1.08 


1.975 


1.919 


0.854 


2.705 


3.318 


-0.81 


7.263 


3.974 


2.962 


7.914 


6.836 


-0.54 


17.05 


7.555 


6.763 


18.54 


12.61 


-0.27 


67.15 


20.68 


25.14 


75.38 


33.84 





179.8 


54.43 


66.41 


155.6 


53.42 


0.27 


69.02 


20.15 


26.44 


56.92 


28.43 


0.54 


16.61 


7.263 


7.305 


16.38 


12.20 


0.81 


5.738 


3.472 


2.526 


5.480 


4.895 


1.08 


2.241 


2.464 


1.031 


2.440 


2.769 


1.35 


1.931 


1.456 


0.887 


1.258 


1.871 


1.62 


1.240 


1.060 


0.384 


0.920 


1.627 


1.89 


0.286 


0.524 


0.0829 


0.607 


1.064 



Table 3. Spatially-measured line fluxes of the prominent emission 
line blends through the 0.35 arcsec wide T-ReCS slit. The detector 
pixels were binned in threes in the spatial direction, giving spaxels 
of 0.27 arcsec. Negative separations are to the South, positive to 
the North. 




10 

Wavelength 



11 12 

i^m) 



inant contribution to the line at 10.88 fim is attributed to 
the 4p-4s transition in Hel, which is not a hydrogenic tran- 
sition. The C in and C IV lines are relatively weak, resulting 
in significant uncertainties which propagate into the con- 
tributions from the higher n carbon lines at 10.5 /xm and 
12.8/xm. About 2/3 of the 12.37 /im and 9.72 /im line fluxes 
are attributable to the He I (7-6) and Hell (10-9) compo- 
nents respectively. [S iv] emission was found to contribute 
only ^15% of the total 10.5-pm emission at the peak posi- 
tion, but increases to ~50% at 0.54" from the centre. The 
recombination line contributions were subtracted from the 
spatially-separated 10.5-/im flux measurements by subtract- 
ing the scaled synthetic spectrum. The resulting spatial pro- 
files are plotted in Fig. [4] 

Because the contribution from [Siv] to the 10.52 /im 
emission blend is small, the peak [S iv] flux is very uncertain, 
This makes the FWHM of the [S iv] emission difficult to esti- 
mate, but it is nonetheless clearly extended with a much flat- 
ter distribution than the continuum emission. [Nell] emis- 
sion was found to contribute ~88±5% of the 12.81-/xm flux 
at the peak, increasing to 95% at 0.54". The emission is 
spatially centrally peaked and symmetrical within the es- 
timated errors, with a FWHM~0.6" and a steeper profile 
than [S iv] . 

3.2 Imaging 

After processing the imaging data, images in the four filters, 
[Siv], [Arm], [Neil] and [Neii]-cont, were produced of 7 Vel 
and a nearby star, SAO 219422, which provided an estimate 
of the PSF. The images were obtained on a number of differ- 
ent nights and in a variety of conditions. The standard star 
images generally showed some low level departures from cir- 
cular symmetry, presumably due to optical aberrations, and 
similar structures were visible in the 7 Vel imaging. The 
PSFs derived from the standard star images have FWHM 
~ 0.4 to 0.6 arcsec. The PSF in both the standard star and 
7 Vel images is typically more centrally concentrated than 
that of the continuum measured from the spectra, which 



Figure 3. The continuum-subtracted spectrum of 7 Vel extracted 
at 0.27" S of the peak (circles connected by the dashed line) 
overlaid by the scaled synthetic spectrum (solid line). The most 
prominent emission line blends are labelled with the wavelength 
and the dominant contributor, but several more components are 
also important in many cases - see Table 2. 



were obtained at higher airmass and in relatively poor see- 
ing conditions; only the images of the WR and standard in 
the [Ar ill] filter show a comparable spatial PSF to the WR 
spectra. 

To obtain images of the spatial distribution of the [S iv] 
and [Nell] emission lines, the underlying continuum and re- 
combination line emission from the WR star had to be sub- 
tracted. This continuum was to have been obtained from 
the [Arm] and [Ne ll]-continuum filter images, which sam- 
ple the continuum emission near to the wavelengths of the 
forbidden-line emission. The [Arm] filter does not include 
any prominent emission lines but the [Ne ll]-continuum fil- 
ter includes a contribution from the Hell 11-10 transition at 
13.12-/im. The contribution from this line to the flux con- 
tained within the [Ne n]-continuum filter was estimated to 
be_~6% using the measurements in lAitken. Roche fc Alien! 
(j 19821 ) and the transmission function of the filter. 

Unfortunately, the continuum subtractions did not ap- 
pear satisfactory, despite using the star spectral energy dis- 
tribution to normalize the continuum fluxes to the wave- 
length of the forbidden emission lines. It seems that the 
variable PSF between the observations dominates and leaves 
substantial residuals, although other factors (e.g. filter leaks 
or ghosts) may also contribute. The spectra indicate that 
while the contribution from the [Nell] emission line to the 
total flux measured within the [Nell] filter is about 15% in 
the central 0.3 arcsec of the image, the contribution from 
the [S iv] line is only about 3%. The approach used in the 
end was to smooth the images and subtract the continuum 
images with the aim of scaling the radial flux to best match 
the spatial flux distribution measured in the spectra. The 
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Position (arcsec) 

Figure 4. Spatially-resolved, extended emission from the [S iv] 
and [No n] forbidden lines in the 7 Velorum outflow obtained from 
the spectra. Unresolved emission from the inner wind is illus- 
trated by the 11/xm continuum and the 12.37/im emission line, 
the latter attributed primarily to the Hel (7-6) line (closed and 
open triangles respectively) and normalised to the [Ne 11] emission 
peak. The extended [S iv] and [Ne n] emission is shown after the 
line deblending described in the text. 



[S iv] subtraction was more troublesome, so the [S iv] image 
of the standard was scaled to the spectrum flux spatial dis- 
tribution and subtracted, but even then the uncertainties 
are considerable. The inner part of the profile is affected by 
the scaling of the continuum image and the outer regions by 
the adopted background level. 

The continuum-subtracted [Ne 11] image, Fig. [5] is cen- 
trally peaked with emission detected to ~ 2 arcsec from the 
centre. The contours are largely centro-symmetric, but with 
slightly enhanced emission to the East. The FWHM of the 
[Nell] image ~0.9", but this figure is rather uncertain be- 
cause of the potential errors in continuum subtraction, and is 
greater than the 0.65 arcsec FWHM obtained from the spec- 
tra. A contour plot of the [S iv] emission is also shown in Fig. 
[5] Extended low- level emission is apparent, again with indi- 
cations of enhanced emission to the East. Note however that 
there appears to be residual emission from the compact core 
in the [Siv] image that is not present in the spatial distri- 
bution of the [S iv] emission derived from the spectra. This 
leads to a compact FWHM comparable to that of the PSF 
standard star at ~0.45 arcsec. Radial profiles generated from 
the [S iv] and [Ne 11] images are displayed together with the 
profile of the comparison star SAO 219422 in the 10.5/im 
filter. Extended emission in the forbidden lines is clearly 
present beyond the core, with the [S iv] profile falling more 
slowly than [Nell] beyond 0.5 arcsec. 

Whilst the detailed forms of the radial profiles are sub- 
ject to substantial uncertainties, all the trial images pro- 
duced retain a high degree of azimuthal symmetry indicating 
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that a spherically symmetric outflow is a good approxima- 
tion. 



4 WIND MODEL 

The spatially extended emission from the fine structure lines 
provides a potential diagnostic of the wind structure. The 
ratio of the [Nell] to [Siv] line intensity probes the wind 
region with densities near the critical densities for the tran- 
sitions, while the radial profiles of these lines will trace the 
level populations, ionic abundances and density of the wind. 

A simple wind model was constructed to compare the 
observed peak and extended emission from [S iv] and [Ne n] 
in the wind with expectations of a freely expanding envelope 
with constant ionization fraction. T wo model winds wer e 
used: (i) a smooth wind as described bv lBarlow et alj l|l988h . 
but modified by the revised distance of 342 pc and a corre- 
sponding mass loss rate of M = 5.7xlO _5 M© yr _1 (rather 
than 460 pc and the mass loss rate of M = 8.8xl0" 5 M o 
yr _1 adopted in iBarlow et al.l (|l988l )): and (ii) a clumped 
wind with a clumping factor / = 10 and a mass loss rate 
of M = 1.5x10 _5 Mq yr _1 , corresponding to the values de- 
termined by de Marco et al (2000) adjusted for the distance 
adopted here of 342 pc. The models calculate the flux admit- 
ted by the T-ReCS slit from a spherical 1/r 2 distribution of 
gas with concentric shells of decreasing electron density, n e , 
away from the core for an optically-thin (at the wavelengths 
of line emission) wind. The instrument slit was projected 
onto the wind model and intervals along the spatial-axis of 
the slit were defined to create bins, representing the pixel 
scale of the observing instrument. The volume of each of 
the individual density shells visible to each bin was then 
calculated. Finally the mean fractional upper atomic level 
population times the radius interval (f u Sr) for each density 
sph ere, shown to trace the observed line flux, from Table 
4 in IBarlow" et al.l l|l988h was multiplied by the normalized 
volume of the corresponding sphere observed in each bin to 
reproduce the flux observed in each bin. 

The collision strengths for [S iv] and [Ne n] adopted are 
8.27 and 0.30 respectively (Saraph & Storey 1999, Johnson 
& Kingston 1987). The Einstein A Values of 7.74 x 10~ 3 and 
8.85xl0 -3 then indicate critical densities n cr it of 3.4xl0 4 
cm -3 for [S iv] and an order of m agnitude greater at 5 xl0 5 
cm -3 for [Nell]. As discussed by Bar low et al.1 0988), line 
formation peaks in the region in the wind where n e ~ n cr i t . 
For the two cases considered here, n cr i t occurs at radii of 
2.7 x 10 15 cm and 7.1 x 10 14 cm for [S iv] and [Ne n] for model 
(i) and at radii of 4.4xl0 15 cm and l.lxlO 15 cm for [Siv] 
and [Nell] for model (ii). 

Because of its lower critical density, the [S iv] emission 
is expected to peak further out in the wind than the [Ne n] 
emission, and for a constant ionization wind the ratio of 
the line fluxes will vary until the density intersected by the 
instrument beam falls below the critical density of [Siv], 
beyond which it should approach a constant value given by 
the atomic parameters and the relative ionic abundances of 
the two species. 

The observed [S iv] and [Ne n] emission line distribu- 
tions given by the models can be compared with the spatial 
distributions of the emission lines extracted from the spec- 
tra. Because they were observed simultaneously, the emis- 
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Figure 5. Top: Contour plots of the [Nell] (left) and [S iv] (right) line intensities. North is up, East to the left. The contours are 
logarithmic in 0.2dex intervals. Bottom: Radial profiles of the comparison star, SAO 219422 at 10.5/xm and the [Siv] and [Nell] line 
intensities in 7 Vel. The solid lines are gaussian fits to the core emission. Because of the uncertainties in subtraction of recombination 
line and continuum emission, the flux levels and distribution in the cores of the line images and profiles may be subject to systematic 
errors 



sion line and continuum distributions extracted from the 
spectrum are more reliable than those obtained from the 
images, which were taken at different times under a range 
of conditions. 

The observed spatial dependence of the [Ne n] to [S iv] 



line intensity ratio is shown together with the ratio calcu- 
lated from the two models in Figure [6] There is reason- 
able qualitative agreement with the observations, although 
the models cannot account for the observed asymmetry and 
neither model matches the observations in detail. In model 
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Figure 6. The [Nell] to [Siv] line ratio (data points) and the line 
ratios calculated from the wind models. The dashed line repre- 
sents model (i)., the unclumped wind, and the solid line represents 
model (ii) the clumped wind with / = 10. 



(ii), the clumping factor of 10 increases the emitting den- 
sity of the wind so that the regions of critical density occur 
at greater distances from the centre, even though the mass 
loss rate is substantially lower. This means that the [Nell] 
to [S iv] line ratio approaches the asymptotic value beyond 
2 arcsec, whereas in the unclumped wind model, the critical 
density is reached at a radius smaller by a factor of 1.6 and 
so the line ratio drops more steeply in the inner region and 
starts to flatten out beyond 1 arcsec. Furthermore, because 
the peak intensity of the [S iv] line is very uncertain, the line 
ratio at the centre has a correspondingly large error, which 
makes the normalization of the fits very uncertain. 

In the outer wind, where the density falls below the 
critical densities of the transitions, the [Nell] to [Siv] line 
intensity ratio approaches 1.0 ±0.2. The dominant ion- 
isation stage of neon is Ne 2+ , so the observations pre- 
sented here cannot place useful constraints on the over- 
all neon abundance. However, they can provide estimates 
of the Ne + to S 3+ ionic ratio. The emission line intensity 
Ia oc Yiv J n u A2,i V, and as the Ne and S lines arise from the 
same volume, the ratio of the line intensities provides an esti- 
mate of the Ne + /S 3+ abundance r atio. The fract i onal u pper 
level populations from Table 3 of iBarlow et al.l (|l98ct ) and 
the older atomic datasets used in that work would have given 
an ionic abundance ratio Ne + /S 3+ ~ 50, while the more re- 
cent atomic datasets for Ne + and S 3+ from the CHIANTI 
database (www.chianti.rl.ac.uk) yield a Ne + /S 3+ ratio of 67 
in the outer wind. 

Fig. [7] compares the spatial distribution of the [Ne n] 
and [S iv] emission lines with the model predictions. The 
spherical wind and density-dependent emission of the for- 
bidden lines does reasonably reproduce the observed emis- 
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Figure 7. Emission observed in the [S iv] and [Ne Ii] lines (data 
points) and the emission distributions calculated from the wind 
models, dashed lines represent the unclumped and solid lines the 
clumped wind profiles. 



sion from a separation of 0.6" outwards, although the range 
of scales sampled is small. The slope of the emission past 
~0.6" appears to confirm that a spherical wind geometry 
with the 1/r 2 density-dependent emission distribution is a 
good match for the observed data in the outer wind. How- 
ever, the models diverge in the inner regions. In the un- 
clumped wind model, the line intensities continue to increase 
steeply towards the central pixel as the critical densities are 
reached relatively close to the centre. The [S iv] emission 
from the clumped wind model starts to flatten at a radius of 
~ larcsec, providing a better match to the observations; the 
model indicates that the [Nell] emission line profile would 
flatten inside a radius of 0.3 arcsec, but unfortunately this 
regime is not well-sampled by the data. 

The clumped wind provides a better overall match to 
the spatial dependences of the emission line intensity distri- 
butions and the [Ne n] to [S iv] line ratio. A modest clump- 
ing factor / ~ 10 provides a reasonable match to the data; a 
more highly clumped wind might improve the fit somewhat, 
but a very highly clumped wind (with / >>10) would flatten 
the profiles further and start to diverge from the observa- 
tions. 

The models made a number of assumptions that could 
have an effect. The emission was ass umed to be o p tically - 
thin at all radii. As discussed by IBarlow et al.l fl988), 
the fine-structure lines are expected to be optically thin 
throughout the outflow, but the infrared free-free continuum 
is expected to become optically thick in the inner dense re- 
gions. The radius at which the continuum becomes optically 
thick can be estimated from the work of Wright and Bar- 
low (1975) who provide an expression for the characteristic 
radius of the emitting region in the outflow. For 7 Vel, the 
characteristic radius at a wavelength of 12 /im for the prop- 
erties adopted for model (i) is r~ 2 x 10 12 cm, and only 
slightly larger for model (ii), so that a negligible fraction 
of the fine-structure line emission will be obscured. Other 
possibilities include a changing ionization structure, possi- 
bly as a result of high ionization in the inner regions, or a 
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non-constant clumping factor. The ionization fraction is not 
expected to vary at such large radii due to the absorption 
of ionizing photons in the inner wind and the decreasing 
recombination rates in the declining-de nsity outer wind re- 
gions ('freeze-out') (|Barlow et alJll988h . but could vary in 
the inner regions. Unfortunately, the data quality presented 
here is not adequate to place stringent limits on these effects, 
but they do indicate that modest degrees of clumping in the 
wind improve the match with observations for the simple 
cases considered here. As shown by Dessert et al (2000), 
adoption of a clumpy wind with lower mass loss rate in- 
creases the ionic abundances of Ne + and S 3+ compared to 
an unclumped wind. The data presented here support the 
need for a clumpy wind. 



5 CONCLUSIONS 

Spatially-resolved imaging and spectroscopy is presented of 
forbidden emission in the outer wind of the Wolf-Rayet star 
in the 7 Velorum binary. The data confirm that forbidden 
[S iv] and [Ne 11] line emission is prominent at large sepa- 
rations (out to ~ 2" radius and ~ 500 x the binary sepa- 
ration) from the star. The images reveal largely symmetric 
outer- wind emission from [Siv] and [Nell], suggesting that 
spherical symmetry is a reasonable approximation for the 
outflow. Spatially-binned spectra are rich in emission lines, 
which clearly follow three distinct spatial emission profiles. 
The He and C lines, and the thermal continuum, all trace 
the distribution of the PSF indicating that the emission is 
unresolved, as expected for recombination lines and free- free 
emission which arise primarily from the dense inner regions 
of the wind. The [S iv] and [Ne 11] emission lines have broader 
spatial profiles. The [S iv] line has a broader spatial emission 
profile than the [Nell] line; this is attributed to its low crit- 
ical density which leads to line formation peaking at larger 
radii in the wind. 

A simple wind model has been constructed, using spher- 
ical shells with a 1/r 2 dependence, to estimate the flux ad- 
mitted by the T-ReCS slit. Two cases were considered: a 
smooth wind with a mass loss rate of M = 5.7xl0 _5 M© 
yr _1 and a clumped wind with clumping factor / = 10 and a 
mass loss rate of M = 1.5x 10 _5 Mq yr _1 . The gradients and 
ratios of the [S iv] and [Ne 11] line profiles at large radii from 
the WR are well-reproduced by the models, showing that a 
spherical wind and 1/r 2 density dependence are appropriate, 
if not unique, solutions. In the inner regions, neither model 
fully reproduces the spatial emission profile of the lines, but 
the unclumped wind model overpredicts the line intensities 
substantially. The clumped wind provides a better fit to the 
data, providing a better match to both the intensity pro- 
files and the spatial distribution of the [Ne 11] to [S iv] line 
ratio. The wind structure could be probed by higher spatial 
resolution observations made under better conditions and 
could be examined further by a similar analysis at mid-IR 
wavelengths of [Nelll] 15.5/mi and [S ill] 18.7 /im using the 
James Webb Space Telescope in the future. 
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